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Abstract 

We study the steady state properties of a phenomenological two-state preda- 
tor model in presence of correlated Gaussian white noise. Based on the cor- 
responding Fokker-Planck equation for probability distribution function the 
steady state solution of the probability distribution function and its extrema 
have been investigated. We show for a typical value of noise correlation there 
is a giant loss of bistability which in turn prevents the prey population from 

going into extinction. 
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I. INTRODUCTION 



The subject of noise-induced transition has got wide apphcations in the field of physics, 
chemistry and biology In most of these theories the noise affects the dynamics through 
system variable, i.e., the noise is multiplicative in nature. The focal theme of these inves- 
tigations is the study of steady state properties of the system where the fluctuations, in 
general, are applied from outside and are independent of the system's characteristic dissipa- 
tion. Such systems are generally termed as open systems , since they lack the principle of 
detailed balance which ensures fluctuation-dissipation relation to hold for the thermodynam- 
ically closed systems. However, it may also happen that the external fluctuations instead 
of affecting only some system's parameters affect the system directly, i.e., they drive sys- 
tem dynamics multiplicatively as well as additively. Because the two noise processes owe a 
common origin they get correlated in the relevant timescale of the problem 0,0. Correlated 
noise processes have found applications in studying steady state properties of a single mode 
laser 0, in analyzing bistable kinetics 0, in giant supression of activation rate 0, in pro- 
ducing directed motion in spatially symmetric periodic potentials in studying stochastic 
resonance in linear systems ||^, in steady state entropy production |jTO|, etc. In this brief 
communication we investigate a simple noise-driven two-state predator model and show 
how noise correlation can dynamically prevent the prey population from extinction. 



II. THE MODEL 

To start with we consider an environment of the prey which in absence of predation grows 
logistically and at the same time its density in a territory depends linearly on a constant 
source of migration. We also consider a population of predators in the given territory which 
lives by feeding on prey. The characteristic time scale over which the population of prey 
and predator varies are very much different, so one can consider the predator population 
to be constant within the generation time of prey. The predators are engaged in two types 
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of activities, viz, hunting or resting. The time scale of predator's two activities are very 
short compared to the generation time of prey, i.e., tr, th ^ fi^^ where tji and th are the 
characteristic average time of resting and hunting, respectively and fi is the birth rate of 
prey. The activity of the predator in the territory resembles the mode of action of enzymes or 
catalysts in a chemical reaction. The enzymes or catalysts in a chemical reaction transform 
substrates in a continuous manner without destroying themselves. The constant predator 
population acts in a similar way by feeding on the prey. To put this ideas in a quantitative 
way we write the evolution equations for the predator and prey 

X = A + fxx[l-^^-^XY , (1) 
Y = - —XY + —Z (2) 

where X is the density of prey in a given territory. The constant A in Eq.(|I|) is due to a 
constant source of prey through immigration. The second term in (|ip is the Fisher logistic 
growth term with birth rate fi and carrying capacity K. Y and Z are the numbers of 
predators in the hunting and resting state, respectively. E is the total constant population 
of the predators, i.e., E = Y(t) + Z{t) = constant. The last term in (|l]) describes the decay 
rate of prey. The model is hybrid in nature in the sense that it has virtue of taking into 
consideration of the logistic growth model as well as of the predator-prey model. 

Following Ref. |jl|] we now consider that the predator population, E is small compared 
to prey population X. To study the overall dynamics within the timescale we make the 
following transformation 

Th = eT*H , = £r* , Y = sY* and Z = eZ* (3) 

where e is a small quantity, t^, rjj are quantities of order fi^^ and Y*, Z* are quantities of 
order X. Using @ in (P and (0) we arrive at 

X = A + ^x(l-|)--ixF* , (4) 
1 1 . 



Now eliminating Y* from (D and using the limit e — we arrive at the following 



dimensionless evolution equation for prey 



X = a + x{l — Ox) — (3 



X 



(6) 



1 + X 



where 



T* 

x=^X 



Art 



E 



and 9 



(7) 



a = 



TrK 



It is interesting to note that the third term in Eq.(P) is the predation term which essentially 
emerges from the two-state of predator activities. The steady state solution of Eq.(^) shows 
a cusp type of catastrophe. The corresponding critical point (ac, Pc, Xc) is given by fl]] 



The necessary condition to have a physically realizable critical point i.e., for Oc, Xc to be 
positive, is 6* < 1. Thus the steady state curve of a; as a function of (3 always shows a bistable 
region for small values of 6. The smallness condition may be maintained by increasing the 
carrying capacity K or by decreasing the ratio th/t^i. 



two quantities which appear in the prey evolution equation as a constant and a multiplicative 
factor, respectively. Expressions for a and (3 in (|^) suggest that they are connected by a 
common parameter /x, the birth rate of the prey. Now if due to some environmental external 
disturbance the birth rate of the prey fluctuates, it is likely to affect both a and (3 in the form 
of additive and multiplicative noises which are connected through a correlation parameter. 
Or in other words the external fluctuations affect the parameter (3 which fluctuates around 
a mean value, thus generating multiplicative noise and at the same time environmental 
fluctuations perturbs the dynamics directly which gives rise to additive noise. As a result 
we have the stochastic differential equation in Stratonovich prescription, 




Eq.(|^) is the starting point of our further analysis. It may be noted that a and (3 are the 




(8) 
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where ^{t) and r]{t) are the stationary Gaussian white noises with the following properties 

m) = m) = , (9) 
(e(t)e(t')) = 2a5(t-t') , (10) 

{r]{t)7]{t')) = 2D6{t - t') and (11) 

mvit')) = mm) = ^x^aDf/Hit - f) (u) 

where A denotes the degree of correlation between noise processes ^(t) and rjit) with < A < 
1. Using the above mentioned noise properties we write down the corresponding Fokker- 
Planck equation (in Stratonovich prescription) for the evolution of probability distribution 
function @,||, 



where 



and 



A(., = a + .(1 - to) - + a-f-- - KoD^j^. (14) 



S(^.*) = « + -7r|^-2A(.Z5)"=^ . (15) 



III. STEADY STATE ANALYSIS AND RESULTS 

Using the zero current condition at the stationary state we derive the stationary proba- 
bility distribution function (SPDF) with and oo as the natural boundaries, 



Psix] = N 



Bix 



exp 



B{x') 



(16) 



where is the normalization constant. Using the explicit forms of A{x) and B{x) we have 
the following explicit forms of SPDF 



Ps{x) = N{1 + x)g'' 2 (x) exp[gix^ + q2x'^ + q^x + qif{x)] 



(17) 



where 



g{x) = a + bx + cx"^ (18) 



f{x) = -2/{b + 2cx) for A = 1 

= {2/y/A) arctan[(6 + 2cx)/VA] for < A < 1 



with 



along with 
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3^ ' 

1-20 be a-P-e + 2 b{l-29) + ae b'^9 



'?2 = ^ — ' ^3 

g4 = a 



2c 2c^ ' c ' 

_ b{2a-P + l) {b^ -2ac){a- (3 -6 + 2) _ b\b^ - 3ac)e 

2^ ^ 2^52 ^ 

a{b^ - 2ac)e - b{b^ - 3ac)(l - 26) 



2c3 



and 



(19) 



a = D , b = 2[D - XiaDf/^] , 

c= D + a-2X{(rDy/^ ,A = AaD{l- X"^) (20) 



_ 2a-/3 + l _ b{a -{39 + 2) jb"^ - ac){l - 29) - ab9 

^ ~ Yc 2^52 ^ 2^3 

6(62 _ ^Q 

■ (21) 

The extrema of SPDF is calculated using the condition A{x) — B'{x) = 0, 

a + x{l- 9x) - - + ][''^^'!' = forO<A<l . (22) 



For zero noise correlation, i.e., for A = the last term of Eq.(P^) vanishes and we have 
the extrema of SPDF for pure multiplicative noise processes For zero correlation the 
additive noise has no extra effects in the steady state dynamics. To illustrate this we have 
plotted extrema of SPDF as a function of (3 in Fig.(l) using the parameters given in [|l|. 
For zero noise correlation the curve shows a sharp minima which decreases on increasing 
A. Similarly, in Fig. (2) we have plotted extrema of SPDF as a function of (3 for different 
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values of additive noise strength D with maximum correlation (A = 1). As the additive 
noise strength increases the well gets flattened and almost vanishes for a large enough value 
oiD. 

In Fig. (3) we show the effect of correlation parameter A on SPDF. For a low value of A 
the SPDF shows the typical bistable region ( see Fig. 3 (a) ) which vanishes for higher values 
of A ( see Fig. 3(b) ). As the value of correlation parameter A increases the peak on the lower 
values of x decreases while for a higher value of A we have a single peak at a higher values 
of X. Since x denotes the prey population, it is clear from Fig. (3) that with the increase 
of A values the prey population recovers from going into extinction. In other words, the 
distribution of prey which was mainly peaked about zero (for a low value of A) signifying 
high extinction rate, moves away from zero with the increase of correlation between noises 
thus favouring the prey's survival. Though Gaussian white noise acting independently and 
multiplicatively favours the extinction of prey the extinction rate decreases drastically 
for a simultaneous perturbation of additive and multiplicative white noise originating from 
a common source, hence connected through a correlation parameter. 

From the expressions of /(x), h and c given in Eqs.(^) and ( PDj ) it is clear that for A = 1.0 
we have always a singular distribution for a = D, since it makes both the parameter b and 
c zero and eventually leads to the divergence of all the g's and u. However this divergence 
can be removed for appreciable difference between the a and D values. In Fig. (4) we 
have plotted the typical behaviour of SPDF for maximum correlation A = 1 which shows 
monotonic decreasing behaviour. In contrast to the behaviour shown in Fig. (3), Fig. (4), 
however, shows the hastening of prey's extinction for a full correlation between additive and 
multiplicative noises. 

In this brief communication we have studied the effect of environmental fluctuation of 
the birth rate of the prey in terms of external correlated noise processes which appreciably 
modify the macroscopic behaviour of a two-state predator model. We have shown how the 
correlation between the two noise processes which owe a common origin may drastically 
prevent the extinction of the prey. 
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FIGURES 

FIG. 1. Plot of extrema of SPDF as a function of (5 for different values of noise correlation A 
using a = 4.5, = 0.1, a = 33.0 and D = 3.0. 

FIG. 2. Same as in Fig.(l) but for different values of additive noise strength D. The other 
parameters are same except A = 1. 

FIG. 3. Plot of Ps{x) against x for different values of noise correlation A using a = 4.5, 9 = 
0.1, (3 = 7.5, a = 3.0 and D = 0.3. (a) For low values of A and (b) for high values of A. 



FIG. 4. Same as in Fig.(3) but for A = 1.0 and D = 2.12. 
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Fig.(1) 




Fig.(2) 
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Fig.(3) 
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Fig.(4) 



